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Daytime Attentuation Rates in the Very Low Frequency 

Band Using Atmospherics* 

W. L. Taylor 

(January 25, 1960) 

Daytime attenuation characteristics have been computed by comparing the amplitude 
spectra of atmospheric waveforms recorded at four widely separated stations. The results 
of these attenuation measurements are presented for the band of frequencies from 3 to 30 
kilocycles per second and involving distances of 1,000 to 10,000 kilometers. It was found 
from these data that attenutation was about 7 to 9 decibels per 1,000 kilometers at 6 kilo- 
cycles per second and decreases to about 1 to 3 decibels per 1,000 kilometers at frequencies 
greater than 10 kilocycles per second. The difference in attenuation rate of west-to-east 
propagation relative to east-to-west propagation was about 3 decibels per 1,000 kilometers 
less for frequencies lower than 8 kilocycles per second and about 1 decibel per 1,000 kilometers 
less from frequencies higher than 10 kilocycles per second. 



1. Introduction 

Considerable electromagnetic energy is radiated 
into the atmosphere from lightning discharges. 
These signals are referred to as atmospherics. The 
spectra of signals from return stroke lightning dis- 
charges attain a peak in the VLF region, usually be- 
tween 5 and 20 kc. The radiation field spectra near 
the source are relatively broad i.e., 10 to 15 kc be- 
tween the frequencies at which the spectral amplitude 
is 3 db below the peak amplitude. As the propaga- 
tion distance increases, the spectra begin to narrow 
such that at distances greater than a few thousand 
kilometers the spectral width is only 3 to 5 kc. The 
spectra of atmospherics simultaneously recorded at 
several locations can be analyzed to determine the 
attenuation rates as a function of frequency and 
direction. 

Many recent experimental and theoretical in- 
vestigations have been conducted for purposes of 
better understanding VLF propagation character- 
istics. 1 The objective of this paper is to present the 
results of a study of daytime attenuation rates at 
VLF using atmospheric waveforms recorded at four 
widely separated locations. 

2. Theoretical Considerations 

The Fourier spectrum of an electromagnetic pulse 
can be represented by the integral 



E(u,d)= ["OM-e-'o'dt 



(1) 



where co is the angular frequency (2irj); d is the 
distance in kilometers; t is the time in seconds, and 
G(t,d) is the instantaneous value of the vertical com- 
ponent of the electric field. Also 



G(t,d)=0 for t<0 and t>r 



(2) 



* Contribution from Central Radio Propagation Laboratory, National Bureau 
of Standards, Boulder, Colo, 
i A selection of representative references is appended. 



where r is the length of the pulse waveform. 
Now (1) can be written 

E(w,d) = \E(^d)\e-^ ( - d) - 



(3) 



The function \E(w,d)\ is the amplitude spectrum and 
the function <t>(u,d) is the phase spectrum as a func- 
tion of frequency at some distance d. The phase of 
the spectrum can be used to determine phase velocity 
of propagation [l]. 2 Only the amplitude value, 
\E(oo,d)\, need be considered for purposes of com- 
puting attenuation presented in this paper. 

Assuming the source is an equivalent dipole, the 
spectrum of the vertical electric field at a great- 
circle distance, d, is given by the mode sum 



\E{u,d)\=A{u,)W, 



(4) 



where A(ui) is the amplitude coefficient dependent 
upon the spectrum of the source, and W is the re- 
sultant of the complex sum of all active waveguide 
type modes [2 to 4]. It is extremely difficult to 
evaluate (4) in a practical situation involving present 
experimental techniques. However, the main con- 
tribution to the spectrum is given by 



\E(o>,d)\~A(a>)[a sin dla]~ 1/2 e~ 



(5) 



where a is the radius of the earth in kilometers and 
a(o)) is the attenuation coefficient in nepers. This 
relation expresses the dominant term in a mode 
series representation of the resultant spectrum, and 
is an equality only when there is one active mode in 
the waveguide. This condition was approximately 
met in the present work as it was in previous work [5]. 
The principal task was to determine the attenua- 
tion as a function of frequency. To accomplish this 
the vertical electric field waveform of the same light- 
rung stroke was recorded at several locations. 
The Fourier spectrum of each waveform was then 



2 Figures in brackets indicate tbe literature references at tbe end of this paper. 
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calculated. Using the subscript i to denote the far 
station and the subscript j to denote the near station, 
such that di^xijj it follows from eq (5) that 



. \E( a),dj)\ f sin (dila) ~\ 
exp [«(„) • (d i -d j )]= m;Ml • L s -nr(^a)J 



(6) 



The attenuation coefficient «'(«) expressed in decibels 
per 1,000 km is given by 

,, \ Tom \E(u,dj)\ 1A1 sind f /a"| 10 3 ,_. 
a'(«) = 2 log ' , r 111 l°g ~ — X/ ' 1 T (7) 

wherea'(co)=a(a>)X8.68X10 3 . 

3. Experimental Procedure 
3.1. Equipment 

The vertical component of the electric field of 
atmospherics was received using a vertical antenna 
and wideband receiving equipment. These atmos- 
pheric waveforms were displayed on a multigun 
cathode-ray tube at two sensitivities and recorded 
on continuously moving photographic film. The 
amplitude response of the waveform channel was 
approximately flat throughout the band-pass and 
sloped to 6-db points at frequencies of 1 and 100 kc. 
The phase response was essentially linear with 
frequency within the band-pass. The cathode-ray 
tube display is proportional to the free-space field 
incident at the antenna. Included in the vertical 
electric field channel was a 23-ft vertical antenna, 
antenna cathode follower, band-pass filters com- 
posed of RC networks, a 24-/xsec delay line, wide- 
band amplifier, and oscilloscope. The overall char- 
acteristics of this channel were adequate to faithfully 
record the VLF components of atmospheric wave- 
forms. 

Calibration was performed by inserting a voltage 
into a dummy antenna driving the cathode follower, 
and adjusting amplifier gain to give the desired 
deflection on the cathode-ray tube. This calibration 
voltage was equal to the expected field strength times 
the antenna effective height, assumed to be one-half 
its physical length. 

Amplitude threshold triggering was used to dis- 
criminate against the numerous low-amplitude at- 
mospherics. Signals that exceeded a selected ampli- 
tude, prior to passing through the delay line, acti- 
vated the sweep circuit of the oscilloscope. All 
signals received with the vertical antenna were 
displayed on the oscilloscope from 24 /xsec before to 
1 msec after the time when the amplitude of the 
atmospheric reached the triggering level. Use of 
the delay line minimized possible loss of the early 
part of the waveforms. 

The direction of arrival of atmospherics was dis- 
played on the same multigun oscilloscope. A pair 
of vertical electrostatically-shielded loop antennas, 
arranged at right angles to each other, loop receivers, 
filters, delay lines, and amplifiers composed the 
direction-finding system. The overall transfer char- 



acteristics of each loop-antenna channel closely 
approximated that of the vertical-antenna channel. 
The differences between the amplitude and phase 
responses of the two loop-antenna waveform channels 
were less than 0.5 db and 6 deg, respectively, over 
the entire band-pass. 

The direction-finding display on the oscilloscope 
was activated for approximately 200 /xsec at the 
beginning of each atmospheric. Thus, the direction 
of arrival was determined from the initial portion of 
the signal, which is essentially vertically polarized. 
Accuracy of determining direction of arrival was 
about ±2° for each station. 

The time of arrival of atmospherics was deter- 
mined to and accuracy of about ± 1 msec. This was 
accomplished by recording on the waveform records 
the standard time signals from stations WWV (or 
WWVH) and timing marks from a secondary fre- 
quency standard which was synchronized with the 
seconds "tick" from WWV. The measured arrival 
times of atmospherics were corrected for propagation 
time of the signals from the standard time station to 
the recording site. 

Four sets of this equipment were operated simul- 
taneously at Boulder, Colo.; Salt Lake City, Utah; 
Palo Alto, Calif.; and Maui, Hawaii, on a number of 
occasions during the summer of 1958. The record- 
ing sensitivities and the triggering amplitudes were 
selected for each station to accommodate atmos- 
pherics originating in the western half of the United 
States. The transfer characteristics of the vertical 
antenna waveform channel were frequently checked 
at each station to insure uniformity of responses 
between stations. It was found that the differences 
between the response characteristics of the equip- 
ment at different stations were less than ±0.5 db 
and ±20°, respectively, in amplitude and phase. 
These differences are small for the purposes of this 
paper. 

Atmospheric waveforms resulting from a particular 
lightning discharge were located on the photographic 
records made of each station by virtue of their times 
of reception and directions of arrival. The locations 
of these discharges were determined by triangulation, 
using the direction of arrival indicated at each 
station. 

3.2. Analysis — East-to-West Propagation 

Waveforms from fifty lightning discharges were ob- 
tained from the data collected on July 2, 1958, 1900 
to 2100 G.m.t.; July 7, 1958, 2300 to 0100 (July 8) 
G.m.t.; and July 9, 1958, 1900 to 2100 G.m.t. At- 
mospherics recorded at each of the four stations from 
twenty of these lightning discharges were analyzed. 
This analysis consisted of (1) scaling the individual 
waveforms recorded at each station to determine the 
functions G(t, d), (see eqs (1) and (2)), (2) computing 
the complex spectra (eq (1)) every kilocycles per 
second from 1 to 100 kc, and (3) computing the 
attenuation (see eq (7)) for each pair of stations, i.e., 
Boulder — Maui, Salt Lake City — Maui, and Palo 
Alto — Maui. The distances from the sources of 
each atmospheric to the near station were between 
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1,000 to 2,400 km, while the distances from the 
sources to Maui (the far station) were between 
4,160 to 5,200 km. Accuracy of determining the 
location of the individual atmospherics was about 
±10 percent of the distance to the near station. 
One set of waveforms, their corresponding ampli- 
tude spectrum, and the computed attenuations 
(a'(w)), representative of those used in this analysis, 
are given in the first four figures. The waveforms, 




400 

time , fi sec 
Figure 1. Representative atmospheric waveforms, east to west. 

shown in figure 1, are reproduced on the same time 
base but different amplitude scales, while the record- 
ing stations and the distances to the source are 
indicated. Note in particular the progression to 
shorter time intervals of the points where each wave- 
form successively crosses the zero amplitude n.xis 
as propagation distance increases. The times in 
microseconds of these quasi half-periods, presented in 
figure 2, show clearly the progression to shorter time 
intervals as propagation distance increases. These 
values, from the representatives set of waveforms, 
agree quite reasonably with those found by other 
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Figure 2. (J nasi half-periods of representative atmospherics, 
east to west. 

workers [6, 7], and the behavior is in accord with 
theoretical investigation of the phenomenon [8]. 
Further elaboration concerning the atmospheric wave 
shapes are beyond the scope of this paper. 

The spectra of these atmospherics are presented in 
figure 3. It should be observed that the frequency of 




TO 100 



FREQUENCY, kc 

Figure 3. Representative spectra, east to 
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the peak amplitude of each spectrum is increasing as 
the propagation distance increases. 

The attenuation in decibel per 1,000 km computed 
from these spectral values are presented in figure 4. 
Values of attenuation for each frequency from each 
set of computations, using Boulder to Maui, Salt 
Lake City to Maui, and Palo Alto to Maui, are in 
reasonably good agreement. Accuracy of determin- 
ing differences in distance between each station of a 
pair and the source was about 7 percent. 




7 10 15 

FREQUENCY, kc 



Figure 4. Representative attenuation, east to west. 

Values of the spectral components of frequencies 
lower than about 4 kc and higher than about 30 kc 
are relatively small, less than 20 db below the peak 
spectral value. Errors in the computed attenuation 
resulting from limitations in the described technique 
may become large as the amplitude of the spectra 
approaches very small values. Therefore, computa- 
tions of attenuation were limited to the VLF region 
of the spectrum (3 to 30 kc). 

The analysis of the waveforms from the other 
lightning discharges give results similar to the 
example presented in figures 1,3, and 4. 

Considerable variability was found between the 
shapes of the attenuation versus frequency curves, 
obtained from the data recorded at each pair of sta- 
tions. Six individual curves of attenuation as a 
function of frequency are presented in figure 5. The 
curves are identified by capital letters (A) and their 
corresponding "zero decibel level" by lower case 
letters (a). The names of the stations and distances 
from the source in kilometers are given with the 
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FREQUENCY, kc 

Figure 5. Individual attenuation, east to west. 

A. Palo Alto, 1,040 km to Maui, 4,160 km. 

B. Salt Lake City, 1,250 km to Maui, 4,480 km. 

C. Salt Lake City, 1,470 km to Maui, 4,700 km. 

D. Boulder, 1,650 km to Maui, 4,300 km. 

E. Palo Alto, 1,840 km to Maui, 4,330 km. 

F. Palo Alto, 2,000 km to Maui, 4,320 km. 



figure. Variations of the computed attenuation 
about some mean value can be seen between 10 to 30 
kc for each of these curves. The small, smooth varia- 
tions of less than ±1 db at frequencies above 10 kc, 
observed in curve B may well be produced by inter- 
ference from higher order modes. Abrupt changes in 
attenuation, below 5 kc in all curves, and observed 
near 12 kc in curve C, 17 kc in curve D and E, and 
29 kc in curve F, are most likely produced by various 
limitations in the recording and analysis techniques. 
These abrupt changes were observed only at those 
frequencies where the spectral amplitude component 
was extremely small for one or both of the wave- 
forms recorded at a pair of stations. Small unwanted 
signals and noise, superimposed on the waveforms, 
and errors in tracing and scaling the waveforms may 
vary the resultant amplitude spectra by a few micro- 
volt-seconds per meter (/^vsec/m) at any frequency. 
Therefore, large errors would be expected at 
those frequencies where the computed amplitude of 
the spectra were small and therefore approaching 
the value introduced through limitations in the 
techniques. 

Errors in the computed attenuation will tend to- 
ward a random distribution. Averaging the attenu- 
ation computed from a relatively large sample of 
independent data should give attenuation values 
that are reasonably accurate. 
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The average attenuation-versus-frequency values 
computed from the data recorded at each pair of 
stations are presented in figure 6. The average loca- 
tion of all lightning discharges used in this study was 
at 51°N to 118°W, i.e., the center of gravity for the 
discharges. The average distance from the lightning 
discharges was 1,220 km to Salt Lake City, 1,470 km 
to Palo Alto, 1,530 km to Boulder. All of these 
propagation paths occurred over land. The average 
distance to Maui was 4,580 km, which was all sea 
water except for 480 km at the source end of the 
average path. It can be seen from these distances 
that the average propagation path, over which these 
attenuation values were computed, was about 82 
percent sea water for the Salt Lake City to Maui 
path, 76 percent sea water for the Palo Alto to Maui 
path, and 74 percent sea water for the Boulder to 
Maui path. The midpoint between this center or 
average location and Maui was approximately 34° N 
to 139° W where the direction of propagation was 
about 230° from geographic north. 8 




7 10 15 

FREQUENCY, kc 

Figure 6. Average attenuation by station, cast to west. 

The average of the three curves is shown in figure 
7. This curve represents experimentally-derived 
attenuation coefficients versus frequency for an ap- 
proximately east-to-west path over predominantly 
sea-water conditions for summer, daytime condi- 
tions. The attenuation values represented by this 
curve may be low compared to an all sea-water path. 
An approximate correction to these values can be 
calculated by assuming the attenuation for a hind 
path is 1 db greater for all frequencies than for a sea- 
water path. Using this approach, it can be shown 
that the attenuation values should be possibly in- 
creased by about 0.3 db/1,000 km. 

3.3. Analysis — West- to-East Propagation 

Attenuation coefficients for this study were com- 
puted from the data recorded at Maui and Palo Alto. 

3 Geomagnetic direction of propagation was about 215° and the dip angle of 
the earth's magnetic field was about 56°. 
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Figure 7. Average daytime alien nation, east to west. 

Only one set of attenuation values were obtained 
from each lightning discharge, in which Maui was 
the close station and Palo Alto was the far station. 
The stations at Salt Lake City and Boulder did not 
record the atmospherics used in this part of the 
analysis because of the high level of atmospheric ac- 
tivity relatively near these stations during each 
observation period. Waveforms from 10 lightning 
discharges OCCliring in the western Pacific Ocean were 
obtained from the data collected on July 2, 1958, 
1900 to 2100 G.m.t., and on July 9, 1958, 1900 to 
2100 G.m.t. Eight of these were selected for the 
nt tenuation study. 

Analysis of these data was conducted in the same 
manner as for the east-to-west attenuation computa- 
tions. One pair of representative waveforms is 
shown in figure 8, reproduced on the same time base. 
Again the times of the zero-axis crossings of the Palo 
Alto waveform, recorded at the more distant station, 
are slightly less than for the Maui waveform. 




Figure 8. 



time , fJL sec 
Representative atmospheric waveforms, west to east. 
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The spectra of these two atmospherics are pre- 
sented in figure 9. The frequency of the peak ampli- 
tude of the Palo Alto spectrum occurs about 500 
cycles higher than the corresponding peak of the Maui 
spectrum. Oscillations of the spectra, occurring be- 
tween about 15 to 40 kc, are probably produced, at 
least in part, by an apparant "phase" shift in the 
waveform between the fifth and sixth half cycles as 
seen in figure 8. Phenomena of this type have been 
reported in the literature [9]. 
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Figure 9. Representative spectra, west to east. 

Attenuation in decibels per 1,000 km computed 
from this pair of atmospherics is shown in figure 10. 
Oscillations in the Palo Alto and Maui spectra, 
mentioned above, show their effect here as oscilla- 
tions in the attenuation curve. 

The average attenuation values, computed from 
the selected pairs of waveforms, are presented in 
figure 11. The accuracy of determining the location 
of individual atmospherics was about ±20 percent 
of the distance to the near station (Maui), but the 
accuracy of determining the differences in distance 
between each station and the source was about ±10 
percent. Approximate average distance from the 
source to Maui was 6,460 km, and from the source to 
Palo Alto was 9,160 Ian. The midpoint between 
the center or average location of the discharges and 
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Figure 10. Representative attenuation, west to east. 
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Figure 11. Average daytime attenuation, west to east. 

Palo Alto was approximately 33° N to 176° W where 
the direction of propagation was about 68° from 
geographic north. 4 This curve represents experi- 
mentally derived attenuation coefficients versus fre- 
quency for an approximately west-to-east path, over 
sea water, for summer, daytime propagation. 

4. Concluding Remarks 

The general shape and the absolute values of the 
east-to-west attenuation-versus-frequency curve in 
figure 7 agrees reasonably well with the experimental 
and theoretical works of others. 5 The attenuation 
values for west-to-east propagation in figure 11 are 
somewhat lower than might be expected for a sharply 
bounded ionospheric model having reasonable values 

4 Geomagnetic direction of propagation was about 58° and the dip angle of 
the earth's magnetic field was about 49°. 

5 A selection of representative references is appended. 
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of reflection height, electron density, and collisional 
frequency. The differences between these two at- 
tenuation curves are shown in figure 12. This repre- 
sents the attenuation for east-to-west propagation 
minus the attenuation for west-to-east propagation. 
A signal propagating from west to east is attenuated 
3 db/1,000 km less between frequencies of about 3 to 
8 kc, and 1 db/1,000 km less between frequencies of 
about 10 to 30 kc, than a signal propagating from 
east to west. The efl'ecl of a decreased attenuation 
at frequencies below about 8 kc would be to enhance 
the smooth oscillatory characteristics of atmospherics 
at great ranges. This characteristic has been ob- 
served by Hepburn [6] and by Chapman and Pierce 
[10] who reported the occurrence of long oscillatory, 
quasi-sinusoidal waveforms from a southwesterly 
direction. 
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Figure 12. Increased alien mil ion for east to west compared 
ivith west to east foi daytime propagation. 

The 1 db/1,000 km difference in attenuation rates 
between 10 and 30 kc is also in agreement with 
various propagation studies using transmissions from 
VLF stations. Attenuation rates [11] computed 
from the data of Round et al. [12], indicate that 
propagation from east to west is characterized by an 
increase in attenuation of 1 or 2 db/1,000 km 
compared with propagation from west to east. Meas- 
urements of Heritage et al. [13], indicates that 
attenuation over the Pacific is slightly higher for 
east- to-west than for west- to-east propagation. Sim- 
ilar results were reported by Crombie [14]. 

It seems reasonable to conclude that propagation 
from west to east is better than from east to west. 
Much work remains to be done if this effect is to be 
understood fully. In particular, attenuation versus 
frequency should be determined as a function of the 
direction of propagation with respect to the earth's 
magnetic field lines. Theoretical investigations also 
should be extended to include those effects. 
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